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ABSTRACT 

The interstellar medium (ISM) seems to have a significant surplus of oxygen which 
was dubbed as the “O crisis”: independent of the adopted interstellar reference abun¬ 
dance, the total number of O atoms depleted from the gas phase far exceeds that tied 
up in solids by as mnch as ~ 160 ppm of 0/H. Recently, it has been hypothesized 
that the missing O could be hidden in micrometer-sized H 2 O ice grains. We examine 
this hypothesis by comparing the infrared (IR) extinction and far-IR emission arising 
from these grains with that observed in the Galactic diffuse ISM. We find that it is 
possible for the diffuse ISM to accommodate 160 ppm of 0/H in ^m-sized H 2 O ice 
grains without violating the observational constraints including the absence of the 
3.1 ftm 0-H absorption feature. More specifically, H 2 O ice grains of radii ~4/rm and 
0/H = 160 ppm are capable of accounting for the observed fiat extinction at ^ 3-8 /im 
and produce no excessive emission in the far-IR. These grains could be present in the 
diffuse ISM through rapid exchange of material between dense molecular clouds where 
they form and diffuse clouds where they are destroyed by photosputtering. 
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1 INTRODUCTION 

The depletion of certain heavy elements from the interstel¬ 
lar gas, known as the “interstellar depletion”, was first no¬ 
ticed by Morton et al. (1973) who found that the gas-phase 
abundances of these elements measured by the Copernicus 
ultraviolet (UV) satellite for interstellar clouds are signifi¬ 
cantly lower than in the Sun. Based on the strong correlation 
between the depletions of these elements and the condensa¬ 
tion temperatures at which they are incorporated into solid 
grains in stellar atmospheres and nebulae. Field (1974) pro¬ 
posed that the elements missing from the gas phase must 
have condensed into dust grains. 

Assuming the interstellar reference abundances to be 
the same as the accepted solar abundances at the time 
(Cameron 1973), Greenberg (1974) compared the observed 
interstellar depletions with the abundances of the dust¬ 
forming elements required to account for the observed vi¬ 
sual extinction. He found that a substantial amount of the 
intermediate-weight elements O, C, and N was unaccounted 
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for by the interstellar gas and dust whereas the heavier ele¬ 
ments Si, Mg and Fe were underabundant. 

Snow & Witt (1996) compiled the elemental abundances 
of early B stars and young F and G stars and found that 
their abundances are just ~ 50-70% of the then accepted 
solar abundances of Anders & Grevesse (1989). They ar¬ 
gued that, due to their young ages, the photospheric abun¬ 
dances of unevolved B stars and young F and G stars are 
more representative of the interstellar composition than the 
4.6-billion-year-old Sun. If the interstellar abundances are 
indeed subsolar, the original question of a surplus of miss¬ 
ing, unaccountable elements posed by Greenberg (1974) then 
turned into a shortage of raw materials for making the dust 
to account for the observed extinction. This was most pro¬ 
nounced for carbon. Whereas published values for the Sun 
range from [G/H]g « 350 ppm to «470 ppm, where ppm 
stands for parts per million, Snow & Witt (1995) derived a 
C/H abundance of 225 ppm for Galactic stars. With the 
gas-phase abundance of [C/HJ^^^ ~ 140 ppm (Cardelli et 
al. 1996) or ~ 100 ppm (Sofia et al. 2011)0 subtracted, the 
remaining C/H abundance of ~ 85 ppm or ~ 125 ppm is in- 


^ Jenkins (2014) argued that the gas-phase C/H abundance of 
Sofia et al. (2011) derived from the strong transition of CII at 
1334 A may be more trustworthy than the previous values mea- 
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sufficient to form the carbonaceous dust required by all dust 
models. This, known as the “C crisis”, still holds when one 
compares the most recent determinations of the solar abun¬ 
dance of C/H«269ppm (Asplund et al. 2009) and early B 
stars of C/H«214ppm (Nieva & Przybilla 2012)0 with the 
latest dust models (e.g., Jones et al. 2013) which all require 
> 200 ppm to be locked up in carbonaceous dust. 
Note that the refined studies of the solar abundances since 
Anders & Grevesse (1989) all led to a gradual, downward re¬ 
vision such that the derived solar abundances were no longer 
substantially higher than that of B stars, and hence the C 
crisis remains unalleviated with the new solar and B-star 
abundances. It is worth noting that Lodders (2003) argued 
that the currently observed solar photospheric abundances 
must be lower than those of the proto-Sun because heavy el¬ 
ements have settled toward the Sun’s interior since the time 
of the Sun’s formation some 4.6 Gyr ago. She further argued 
that the protosolar abundances are more representative of 
the solar system elemental abundances. 

The depletion of oxygen is also problematic. Unlike G 
of which the depletion is insufficient to account for the ob¬ 
served extinction, the major solid-phase reservoirs of O in 
the diffuse interstellar medium (ISM) — silicates and metal 
oxides — are insufficient to account for the total O/H abun¬ 
dance missing from the gas phase. Jenkins (2009) examined 
the relative proportions of 17 individual elements that are 
incorporated into dust along 243 different Galactic sight¬ 
lines. He found that the depletion of O in the diffuse ISM 
far exceeds the consumption of O by silicates and metallic 
oxides, except for low-density regions with little depletions. 
This is insensitive to the adopted reference abundance since 
the approach taken by Jenkins (2009) was based on differen¬ 
tial depletions rather than absolute values. Whittet (2010b) 
investigated the depletions of O over a wide range of envi¬ 
ronments from the tenuous intercloud medium and diffuse 
clouds to dense clouds where H 2 O ice is present. He found 
that as much as ~ 160 ppm of O/H is unaccounted for at the 
interface between diffuse and dense phases, again, indepen¬ 
dent of the choice of reference abundances. This surplus of 
O, dubbed as the interstellar “O crisis” by Whittet (2010a), 
poses a severe challenge to our understanding of interstellar 
dust. 

The nondetection of the 3.1/rm 0-H stretching ab¬ 
sorption feature of H 2 O ice in the diffuse ISM rules out 
submicrometer-sized H 2 O ice grains as a significant reser¬ 
voir of O (Whittet et al. 1997, Poteet et al. 2015). Jenkins 
(2009) suggested that large amounts of O could be hidden 
in H 2 O ice grains (or large grains with thick mantles of H 2 O 
ice) that have diameters of the order of or greater than 1 /rm. 
Grains this large will have the 3.1 ^m 0-H feature substan¬ 
tially suppressed^ More recently, Poteet et al. (2015) an¬ 
alyzed the ~ 2.4-36 /rm absorption spectrum of the sight¬ 
line toward Ophiuchi obtained with the Short Wavelength 

sured from the weak intersystem absorption transition of CII] at 
2325 A. 

^ Poteet et al. (2015) found that the B-star Si and Mg abundances 
are not enough to account for the 9.7 grci Si-O absorption feature 
observed toward Ophiuchi (see their Figure 8). 

® As illustrated in Figure 4 of Poteet et al. (2015), for 2.0-3.2 pm- 
sized H 2 O ice grains, with increasing grain size, both the scatter¬ 
ing and absorption at the 3.1 pm O-H stretch decrease. 


Spectrometer (SWS) on board the Infrared Space Observa¬ 
tory (ISO) and the Infrared Spectrograph (IRS) on board 
the Spitzer Space Telescope. They determined the elemental 
abundances of O, Mg, and Si in silicates. Along with the 
upper limit estimates of O in other materials, they found 
that as much as ~ 156 ppm of O/H is unaccounted for along 
the line of sight toward ^ Ophiuchi, a prototypical cool dif¬ 
fuse cloud and an environment near the diffuse-dense ISM 
transition. They argued that the missing reservoir of O must 
reside on very large, micrometer-sized grains (e.g., with radii 
aA3.2 pm) which are nearly opaque to infrared (IR) radia¬ 
tion. 

The aim of this work is to test the hypothesis of pm- 
sized H 2 O ice grains as the repository of the unaccounted 
for ~ 160 ppm of O/H in the diffuse ISM. Roughly speak¬ 
ing, dust scatters and absorbs starlight most effectively 
at a wavelength A comparable to its spherical radius a: 
4a(n — 1)/A ~ 1, where n(A) is the real part of the index 
of refraction of the dust. This can be understood in terms 
of elementary optics: the maxima of the extinction could 
be considered as being caused by the destructive interfer¬ 
ence between the incident and forward-scattered light. The 
phase difference between a ray that traverses a large trans¬ 
parent sphere without deviation (i.e., the forward-scattered 
ray) and a ray that traverses the same physical path outside 
the sphere is Acf = (27r/A)2a(n — 1). The maxima of the 
extinction occur at A(j> = (2j -|- l)7r, where j is an integer, 
with j — 0 corresponding to the principal extinction peak 
(see Bohren & Huffman 1983). 

For H 2 O ice, n ~ 1.3 at AAlpm (Warren 1984), thus, 
pm-sized H 2 O ice grains will cause extinction at A > 1 pm. 
Therefore, we will examine in ^Uthe consistency between the 
observed IR extinction and the extinction resulting from the 
standard silicate-graphite model with the addition of a pop¬ 
ulation of pm-sized H 2 O ice grains of O/H = 160 ppm. On 
the other hand, H 2 O ice grains are transparent in the ul¬ 
traviolet (UV) and visible wavelength range and hence they 
will be cold and emit in the far-IR at AA;200 pm. Therefore, 
we will examine in 321 the IR radiation emitted by pm-sized 
H 2 O ice grains and make a comparison between the model 
IR emission and that observed for the Galactic diffuse ISM. 
Their origin and survival will be discussed in fJH The major 
conclusions are summarized in SjS] 


2 EXTINCTION 

For H 2 O ice grains to be considered as a viable component 
of the diffuse ISM, they should satisfy the constraints placed 
by the observed extinction from the far-UV to the IR, in¬ 
cluding the absence of the 3.1 pm absorption feature of H 2 O 
ice. The observed UV/optical interstellar extinction can be 
characterized by a single parameter Rv (Gardelli et al. 1989, 
CGM)0 As elaborated in Wang, Li & Jiang (2014), while the 
UV/optical extinction can be closely fitted by the classical 
silicate-graphite model (Mathis et al. 1977, Draine & Lee 
1984), this model predicts a power-law of Ax oc A“^'^® at 

^ Rv = Ayf (Ag — Ay) is the total-to-selective extinction ra¬ 
tio, where Ag and Ay are the blue- and visual-band extinction 
at A^ = 4400 A and Ay = 5500 A, respectively. For the Galactic 
average, Ry fa 3.1. 
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Figure 1. Fitting the Ry = 3.1 extinction curve from the UV/optical to the near- and mid-IR with (1) amorphous silicate (green short 
dashed line), (2) graphite and PAHs (blue dashed line), and (3) /rm-sized H 2 O ice grains (magenta dash-dot-dotted line) of O/H = 160 ppm 
with a log-normal size distribution characterized by ao ~ 4.0 /rm and a fs 0.4. The thick red solid line plots the model fit which is the 
sum of silicate, graphite/PAHs, and /rm-sized H 2 O ice. The black solid line plots the sum of silicate and graphite/PAHs. The symbols 
plot the observed extinction: the cyan open triangles plot the i?v = 3.1 UV/optical/near-IR extinction, and the other symbols plot the 
mid-IR extinction (see text). 


1 fj,m < A < 7 fim (Draine 1989) which is too steep to be 
consistent with the subsequent observations made by ISO 
and Spitzer since numerous observations suggest that the 
mid-IR extinction at 3 pm < A < 8 pm is flat or “gray” for 
both diffuse and dense environments (see Figure [TJ, includ¬ 
ing the Galactic center (Lutz 1999, Nishiyama et al. 2009), 
the Galactic plane (Indebetouw et al. 2005, Jiang et al. 2006, 
Gao et al. 2009), the Goalsack nebula (Wang et al. 2013), 
and nearby star-forming regions (Flaherty et al. 2007). All 
these observations appear to suggest an “universally” flat 
mid-IR extinction law, with little dependence on environ¬ 
ments. 

We construct an “observed” extinction curve for the 
diffuse ISM as follows: (i) for 0.125/rm < A < 3 pm, we 
take the Galactic average of Rv = 3.1 as parameterized by 
CGM and divide it into 95 wavelengths, equally-spaced in 
In A; (ii) for the mid-IR extinction at 3 pm < A < 8 pm, we 
first obtain a weighted “average” from the observed extinc¬ 
tion shown in Figure [U with twice as much weight given to 
the diffuse sightlines toward the Galactic center^ we then 

® The extinction along the line of sight toward the Galactic cen- 


interpolate the “average” mid-IR extinction into 25 logarith¬ 
mically equally-spaced wavelengths. 

We aim at reproducing the observed extinction from the 
UV/optical to the near- and mid-IR with a mixture of amor¬ 
phous silicate dust and carbonaceous dust as well as H 2 O ice 
grains. The UV/optical extinction is predominantly caused 
by sub-^m- and nano-sized grains (Li 2004). The nonde- 

ter is believed to be dominated by dust in the diffuse ISM as 
revealed by the detection of the 3.4 pm aliphatic C—H absorption 
feature which is absent in dense regions (Sandford et al. 1991, 
Pendleton et al. 1994, Tielens et al. 1996). However, it does con¬ 
tain molecular cloud materials as revealed by the detection of the 
3.1 and 6.0 pm H 2 O ice absorption features, e.g., the sightline to¬ 
ward the Galactic center source Sgr A* suffers ~ 30 mag of visual 
extinction (McFadzean et al. 1989), to which molecular clouds 
may contribute as much as ~ 10 mag (Whittet et al. 1997). As 
shown in Figure [T] the ~3—8 pm extinction curve of the Galactic 
center (Lutz 1999, Nishiyama et al. 2009) is less flatter than that 
of other regions. When constructing the “observed” extinction 
curve, if one assigns less weight to the Galactic center sightlines, 
one would obtain an even flatter ~ 3—8 pm extinction law. There¬ 
fore, even more O/H could be tied up in pm-sized H 2 O ice grains. 
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tection of the 3.1 /im absorption feature of H 2 O ice in the 
diffuse ISM places an upper limit of 2 ppm (Whittet et al. 
1997) and ~9ppm (Poteet et al. 2015) of O/H in sub-/rm- 
sized H 2 O ice grains toward Cyg OB2 No. 12 and C) Oph, 
respectively. Therefore, the observed UV/optical extinction 
is mainly produced by the silicate and carbonaceous compo¬ 
nents while the H 2 O ice component, if present in the diffuse 
ISM, must be larger than ~ 1 /rm and cause extinction at 
A.^1 pm. 

Following Weingartner & Draine (2001; WDOl), we 
model the observed extinction in the wavelength range of 
0.125 pm < A < 8 pm in terms of the silicate-graphite- 
PAH model combined with a population of pm-sized H 2 O 
ice grains. We take the size distribution functional form of 
WDOl for the silicate and graphitic component, assuming 
the latter extends from grains with graphitic properties at 
radii a^50A, down to grains with PAH-like properties at 
very small sizes (see Li & Draine 2001). The dielectric func¬ 
tions of amorphous silicate, graphite, and H 2 O ice are taken 
from Draine & Lee (1984) and Warren (1984). The WDOl 
model employs two log-normal size distributions for two 
populations of PAHs which respectively peak at flop, 00,2 
and have a width of cri, (J 2 , consuming a C abundance of 
5c. 1 , 5c ,2 (per H nuclei). Following Draine & Li (2007), we 
adopt ao.i = 3.5 A, cti = 0.40, 5c ,1 = 45 ppm, 00,2 = 20 A, 
(72 = 0.55, and 5c ,2 = 15 ppm. These parameters were con¬ 
strained by the observed near- and mid-IR emission. For the 
pm-sized H 2 O ice component, we also adopt a log-normal 
size distribution of peak size no and width cr: 


1 dn 
riH da 


"(27r)3/2 
1 

X - exp < — - 
a ' ^ 


^ exp ( 4.5cr ) ^ 5icEPicEmH 

picEfflo'^ 2 


ln(a/ao) 


( 1 ) 


where mn is the atomic H mass, picE and picE are respec¬ 
tively the mass density and the molecular weight of H 2 O ice 
(picE ~ 1.0gcm“3 and picE « 18), and 5ice is the O abun¬ 
dance per H nuclei locked up in pm-sized H 2 O ice grains. 

Following WDOl, we use the Levenberg-Marquardt 
method (Press et al. 1992) to minimize the fitting error 
between the observed and model extinction (see Wang, 
Li & Jiang 2015). As shown in Figure [T] together with 
silicate, graphite and PAHs, H 2 O ice grains with ao « 
(4.0 ± 1.0) pm, a « (0.4 ± 0.1), and 5ice = 160ppm of O/H 
satisfactorily reproduce the ~ 3-8 pm mid-IR extinction. 
Figure [T] also shows that these pm-sized H 2 O ice grains are 
“gray” in the UV/optical and make negligible contribution 
to the observed UV/optical extinction. The size distribu¬ 
tions of the silicate and graphite components which closely 
reproduce the observed UV/optical extinction are character¬ 
ized by the following parameters: Cg « 5.75 x 10~^^, Og « 
— 1.40, /3g « 0.0291, at,g « 0.00818 pm, and Oc.g « 0.173 pm 
for graphite, Cs « 7.56 x 10“'^'^, as « —2.19, /3s ~ —0.586, 
fflt.s ~ 0.204 pm, and Oc.s = 0.1pm for silicate (see WDOl 
for the definition of each parameter). 

The model extinction does not show the 3.1 pm absorp¬ 
tion feature which is absent in the diffuse ISM. However, 
it exhibits a narrow, minor structure at ~ 2.8 pm arising 
from the scattering of the O-H stretch of H 2 O ice. The ex¬ 
tinction toward the Galactic center derived by Lutz (1999) 


based on the recombination lines of atomic H was not suffi¬ 
ciently resolved in wavelength to rule out this structure. We 
note that Whittet et al. (1997) reported the detection of a 
shallow feature centered at 2.75 pm in the near-IR spec¬ 
trum of Cygnus OB2 No. 12 obtained by the Kuiper Air¬ 
borne Observatory (KAO) and ISO/SWS. They tentatively 
attributed this feature to the OH groups of hydrated sili¬ 
cates. However, it was later found to be an artifact caused 
by the calibration uncertainty (see Whittet et al. 2001). 


3 INFRARED EMISSION 

To examine whether the model, with the inclusion of a pop¬ 
ulation of pm-sized H 2 O ice grains, will result in too much 
emission in the far-IR, we calculate the equilibrium tem¬ 
peratures of silicate and graphite grains of radii larger than 
~250A as well as pm-sized H 2 O ice grains heated by the 
Mathis, Mezger, & Panagia (1983, MMP83) interstellar radi¬ 
ation field (ISRF; see Figure [21l. We also calculate the tem¬ 
perature probability distribution functions of PAHs, small 
graphite and silicate grains of radii smaller than ~ 250 A 
since the heat contents of these ultrasmall grains are smaller 
than or comparable to the energy of a single UV photon and 
therefore, they will be transiently heated and will not attain 
an equilibrium temperature (see Draine & Li 2001). 

For a given size, H 2 O ice grains are colder than silicate 
and graphite grains (e.g., T « 14.3, 15.5, 19.7K for H 2 O ice, 
silicate and graphite grains of a = 0.1pm, respectively). 
For a given composition, the dust temperature decreases as 
the grain size increases (see Figure [2}. The best-fit H 2 O ice 
grains of radii of ~4pm have an equilibrium temperature 
of ~8.9K. As shown in Figure O these grains mainly emit 
at A^200 pm. 

Figure |3] presents the model IR emission which com¬ 
bines the contributions from silicate grains, graphite/PAH 
grains, and H 2 O ice grains. We compare the model IR emis¬ 
sion with the diffuse ISM observed by Planck and the Dif¬ 
fuse Infrared Background Experiment (DIRBE) and the Far 
Infrared Absolute Spectrophotometer (FIRAS) instruments 
aboard the Cosmic Background Explorer (COBE). Figure [3] 
shows that the model closely fits the observed emission from 
the near-IR to the far-IR and millimeter (mm). The pm- 
sized H 2 O ice component is “gray” in the UV/optical (see 
Figure [TJ and does not absorb much and therefore by im¬ 
plication, it does not emit much in the far-IR. Although 
the mass of the pm-sized H 2 O ice component exceeds that 
of graphite/PAHs by ~7%, the UV/optical extinction con¬ 
tributed by the former is smaller than that of the latter 
by a factor of ~50. Moreover, the UV/optical extinction of 
H 2 O ice grains is dominated by scattering, instead of ab¬ 
sorption. The silicate-graphite/PAH-ice model results in a 
total IR intensity of ~4.48 x 10“^'* ergs”*^ H“^. The frac¬ 
tional contributions of silicate, graphite/PAHs, and H 2 O ice 
are approximately 27.9%, 71.7%, and 0.34%, respectively. 


4 DISCUSSION 

We have shown in ^ and ^ that the diffuse ISM has no 
problem in accommodating the unaccounted for ~ 160 ppm 
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Figure 2. Equilibrium temperatures for graphite (red solid line), silicate (blue dashed line), and H 2 O ice grains (orange dash-dotted 
line) heated by the MMP83 ISRF. 


of 0/H. With the missing ~ 160 ppm of 0/H tied up in fim- 
sized H 2 O ice grains, the silicate-graphite/PAH-ice model 
closely reproduces the observed extinction from the far-UV 
to the mid-IR (see Figure [1]) as well as the observed ther¬ 
mal emission from the near-IR to the submm/mm wave¬ 
lengths (see Figure El). However, the “C crisis” still per¬ 
sists. The model requires « 225 ppm to be de¬ 

pleted in graphite/PAHs (with ~ 60 ppm of C/H in PAHs, 
see Li & Draine 2001). With the gas-phase C/H abundance 
of [C/H]^^,, Ks 140 ppm (Cardelli et al. 1996) or [C/H]g^,, « 
100 ppm (Sofia et al. 2011) included, the total required C/H 
abundance becomes C/H « 365 ppm or C/H « 325 ppm, ex¬ 
ceeding that of solar ([C/H]g « 269±31 ppm, Asplund et al. 
2009), proto-Sun ([C/H]q « 288 ± 27 ppm, Lodders 2003), 
and early B stars ([C/H]^'« 214±20ppm, Nieva& Przybilla 
2012). We note that the shortage of C/H is true for all dust 
models (e.g., Zubko et al. [2004; ZDA] required [C/H]^^^^ « 
244ppm, Jones et al. [2013] required [C/H]^,^^^ « 233ppm). 
Mathis (1996) argued that porous dust is more effective in 
absorbing starlight than compact dust and therefore one 
might require fewer C atoms to account for the observed ex¬ 
tinction. However, Li (2005a) showed that, through an anal¬ 
ysis based on the Kramers-Kronig relation of Purcell (1969), 
models invoking porous dust cannot appreciably lower the 
C/H consumption. A likely solution to the C crisis prob¬ 
lem is that the stellar photospheric abundances may be con¬ 


siderably lower than that of the interstellar material from 
which young stars are formed (see Li 2005a). Alternatively, 
if the interstellar C/H abundance is like the current proto¬ 
solar abundance of C/H~339ppm (Asplund 2009) when 
taking into account the Galactic chemical enrichment over 
the past 4.6Gyr (Chiappini, Romano, & Matteucci 2003), 
with [C/Hjg^,, ~ 100 ppm (Sofia et al. 2011) subtracted, the 
available C/H abundance (~ 239 ppm) is approximately con¬ 
sistent with the latest dust models and the so-called C-crisis 
will be alleviated. 

The silicate-graphite/PAH-ice model requires 
[Si/H] 

dust ^ [Mg/H],,,,, « 40.4 ppm to be depleted in 
silicate dust. Here we assume a stoichiometric composition 
of MgFeSi 04 for the silicate dust which depletes all the 
Si, Mg, and Fe atoms. For Si, it is higher than that of 
solar ([Si/H]g « 32.4±2.2ppm, Asplund et al. 2009) and 
early B stars ([Si/H]^ « 31.6 ± 3.6ppm, Nieva & Przybilla 
2012), it is consistent with the proto-Sun abundance of 
Si/H « 40.7±1.9 ppm (Lodders 2003). In contrast, the Jones 
et al. (2013) model consumed [Si/H]^^^^ « 50 ppm. While 
the ZDA model used fewer Si atoms ([Si/H]^,^^^ « 36 ppm) 
than the WDOl model ([Si/H]j^,jj « 48ppm), the former 
produces much elss extinction at A > 1 /rm than the latter 
(see Figure 23.11 in Draine 2011). 

Historically, ices were among the first grain species sug¬ 
gested to be present in the ISM and responsible for the 
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Figure 3. Comparison of the observed emission from the diffuse ISM (black diamonds, black squares, and orange circles) to the 
model which is the sum of silicate (green short dashed line), graphite/PAHs (blue dashed line), and /rm-sized H 2 O ice grains (magenta 
dash-dot-dotted line). Red triangles show the model spectrum (red solid line) convolved with the DIRBE filters. Observational data are 
from DIRBE (black diamonds; Arendt et al. 1998), EIRAS (black squares; Finkbeiner et al. 1999), and Planck (orange circles; Planck 
Collaboration XVII 2014). 


observed extinction (see Li 2005b). Oort & van de Hulst 
(1946) proposed the “dirty ice” model consisting of satu¬ 
rated molecules such as H 2 O, CH 4 , and NH 3 . In the 1960s, 
the first attempt to search for the 3.1pm feature of H 2 O 
ice in the diffnse ISM was unsuccessful (Danielson et al. 
1965, Knacke et al. 1969). This led to the abandonment of 
the “dirty ice” model and the proposition of “organic re¬ 
fractory” resulting from the photoprocessing of “dirty ices” 
as an interstellar grain component (Greenberg et al. 1972, 
1995)0 Based on the then “surplus” of O, C, N elements, 
Greenberg (1974) speculated that these “excessive” elements 
could be bound in “interstellar snowballs” ranging in size 
from baseballs to comets. To the best of our knowledge, the 
ZDA model is the only contemporary model in which sub- 


pm-sized H 2 O ice was included as a grain component and 
accounted for ~4.4% of the total dust mass. 

A crucial issue for the hypothesis of pm-sized H 2 O ice 
grains as the reservoir of the missing O atoms is how they 
form and survive in the ISM. In dense molecular clouds, H 2 O 
forms through hydrogenation of O on the surfaces of sub- 
pm-sized silicate and carbonaceous grains and in this waw 
a thick ice mantle could be built up (Boogert et al. 2015)0 
In contrast, ice mantles cannot be built up in the diffuse 
ISM since they will be rapidly removed by photosputtering 
(Barlow 1978a,b; Draine & Salpeter 1979). We note that 
in the diffuse ISM, not only H 2 O ice but also silicate and 
graphite are destroyed at a rate faster than their production 
(McKee 1989). This led Draine (1990) to conclude that the 


® The core-mantle model of Li & Greenberg (1997) assumes that 
the silicate core is coated by a mantle of organic refractory mate¬ 
rial. Whittet (2010b) suggested that the organic refractory matter 
may potentially account for some of the missing O atoms. 


^ One would expect an ice grain to have a sub-pm-sized nu- 
cleation core of silicate or graphite. But for the pm-sized ice 
grains of interest here, the effects of the core on the extinc¬ 
tion and emission are negligible since the core only accounts for 
~ (0.1/4)® ^ 1.6 X 10“® of the volume of a grain of a = 4 pm. 
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bulk of the solid material in grains actually condensed in 
the ISM rather than in stellar outflows. Draine (1990) ar¬ 
gued that there must be rapid exchange of matter between 
the diffuse ISM and molecular clouds since the bulk of grain 
growth can proceed rapidly only in dense regions. Depend¬ 
ing on the posited mass exchange scenario, this implies a 
turnover time (tmc) of 3 x 10®-2 x 10^ yr for molecular 
clouds (Draine 1990). If the photosputtering lifetimes (rpd) 
of /rm-sized H 2 O ice grains are not shorter than tmc, the ice 
grain model might be viable: before they are destroyed by 
photosputtering, very large ice grains in the diffuse ISM are 
continuously replenished by the freshly-condensed ice grains 
formed in dense molecular clouds through accretion and co¬ 
agulation. 

Let ypd(A) be the photosputtering yield of H 2 O ice (i.e., 
the number of H 2 O molecules desorbed per absorbed UV 
photon of wavelength A). We derive the photodesorption rate 
of H 2 O ice as following: 

N = dN/dt = [ d\ , (2) 

J hc/X 


where h is the Planck constant, c is the speed of light, 
Cabs (a. A) is the absorption cross section of H 2 O ice grain 
of radius a at wavelength A, and Jx is the intensity of the 
MMP83 ISRF. We calculate the photosputtering lifetimes of 
H 2 O ices of radii a from 


Tpd(a) 


1 AttC piCE 
N 3picEmH 


(3) 


Westley et al. (1995) experimentally measured the pho¬ 
todesorption yield of H 2 O ice exposed to Lyman-a pho¬ 
tons (10.2 eV or 1216 A) at temperature T to be Ypd ~ 
Yq -b Yi exp (—i?/fcr), where k is the Boltzmann constant, 
Yo = 0.035 ± 0.002, Yi = 0.13 ± 0.10, and 15 = (29 ± 6) x 
10“^ eV. For H 2 O ice grains of radii a = 4 pm and temper¬ 
ature T = 8.9 K (see Figure [2J, we estimate the photosput¬ 
tering lifetime to be Tpd ~ 5.8 x 10® yr in the diffuse ISM. 
This lifetime may be underestimated since the photodesorp¬ 
tion yield of Westley et al. (1995) could have been over¬ 
estimated as their data were obtained through laser beam 
irradiation which might have induced local point heating 
and thus adding some sublimation to the pure sputtering ef¬ 
fect. More recently, Oberg et al. (2009) experimentally deter¬ 
mined the photodesorption yields of H 2 O ice using a hydro¬ 
gen discharge lamp (~ 7-10.5 eV). They derived an average 
photodesorption yield of Vpd ~ 1-3 x 10“® + 3.2 x 10“® x T. 
This gives a photosputtering lifetime of Tpd ~ 7.4 x 10® yr for 
H 2 O ice grains of a = 4 pm at T = 8.9 K. Moreover, Anders- 
son et al. (2006) modeled the photodesorption process at the 
atomic level by simulating the interaction between photons 
of ~ 8-9.5 eV and ice surfaces. They estimated the photodes¬ 
orption yield to be Tpd « 4 x 10“"^ for amorphous ice. This 
leads to a photosputtering lifetime of Tpd ~ 2.9 x 10^ yr for 
H 2 O ice grains of a = 4 pm. Therefore, it seems likely that 
the diffuse-phase pm-sized H 2 O ice grains are continuously 
replenished on a timely fashion by the dense-phase materials 
through the turnover of molecular clouds within a timescale 
of TMC ~ 3 X 10®-2 X 10^ yr before they are destroyed by pho¬ 
tosputtering in the diffuse ISM. In contrast, a sub-pm-sized 
ice grain (coated on a silicate or graphite core) responsible 
for the 3.1pm absorption feature in dense clouds will be 
quickly removed by photosputtering in diffuse clouds at a 


rate faster than that of an a = 4 pm ice grain by a factor of 
-(4/0.1)® « 6.4 X 10'*. 

Finally, we note that the “missing” O atoms do not all 
have to be tied up in H 2 O ice grains. As mentioned earlier, 
some of the “missing” O atoms could be hidden in organic re¬ 
fractories (Whittet 2010b). On the other hand, the observed 
flat extinction at — 3-8 pm could also partly be caused by 
other dust components, e.g., pm-sized graphite (see Wang, 
Li, & Jiang 2015). 


5 SUMMARY 

While the ISM seems to be short of the element C to form 
a sufficient amount of carbonaceous dust (together with sil¬ 
icate dust) to account for the observed extinction, the ele¬ 
ment O seems overabundant and as many as — 160 O atoms 
(per 10® H atoms) are unaccounted for by their presence in 
gas and dust. We have examined the hypothesis of pm-sized 
H 2 O ice grains as the reservoir of the missing O atoms. It 
is found that the diffuse ISM has no difficulty in accommo¬ 
dating — 160 ppm of O/H in pm-sized H 2 O ice grains, con¬ 
firming the earlier suggestions made by Jenkins (2009) and 
Poteet et al. (2015). With a radius of —4 pm, these grains 
are “gray” in the UV/optical and contribute very little to the 
UV/optical extinction. The 3.1pm 0-H stretching feature 
of these grains is significantly suppressed, consistent with 
the nondetection of this feature in the diffuse ISM. They 
absorb and scatter effectively in the mid-IR and are capable 
of accounting for the observed flat extinction at — 3-8 pm. 
Being relatively transparent in the UV/optical, they do not 
absorb much and therefore they are cold and mainly emit 
in the far-IR at AA,200pm, accounting for only —0.34% of 
the total IR power of the Galactic diffuse ISM. With a pho¬ 
tosputtering lifetime of Tpd — 5.8 x 10®-2.9 x 10^ yr longer 
than or comparable to the turnover timescale of molecular 
clouds of Tmc — 3 x 10®-2 x 10^ yr implied by the observed 
large depletions of Si and Fe elements in the diffuse ISM, 
H 2 O ice grains of radii aA,4pm could be present in the dif¬ 
fuse ISM through rapid exchange of material between dense 
molecular clouds where they form and diffuse clouds where 
they are destroyed by photosputtering. 
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